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Abstract 

Simple random walk considerations are used to interpret rodent population data 
collected in Hantavirus-related investigations in Panama regarding the short-tailed 
cane mouse, Zygodontomys brevicauda. The diffusion constant of mice is evaluated 
to be of the order of (and larger than) 200 meters squared per day. The investigation 
also shows that the rodent mean square displacement saturates in time, indicating 
the existence of a spatial scale which could, in principle, be the home range of the 
rodents. This home range is concluded to be of the order of 70 meters. Theoretical 
analysis is provided for interpreting animal movement data in terms of an interplay 
of the home ranges, the diffusion constant, and the size of the grid used to monitor 
the movement. The study gives impetus to a substantial modification of existing 
theory of the spread of the Hantavirus epidemic which has been based on simple 
diffusive motion of the rodents, and additionally emphasizes the importance for 
developing more accurate techniques for the measurement of rodent movement. 
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1 Introduction 



The Hantavir us epidemic i s of great concern to hum a n health in many region s 
of the world ([Yates et all 120021: iMills et all Il99flt IParmenter et all Il99flh . 
The discovery of Hantavirus in the New World took place after an outbreak 
of a severe disease in the region of the Four Corners in the North American 
Southwest, in 1993. The agent, Hantavirus Sin Nombre, is carri ed mainly by 



the extremely common d eer mouse, Peromyscus maniculatus (jNichol et al. 



199.4 IChilds et all Il994h . Since the discovery an enormous effort has been 
devoted to understanding the ecology and epidemiology of the virus-mouse 



assoc i ation, with the ul t imate goal being prediction of human risk ((Mills et al. 



1999; IParmenter et all Il998f ). Numerous species of the virus are known in 
the Americas, ea ch one of them almost e xclusively associated with a single 
rodent reservoir ([Schmaliohn et all [1997). Human disease caused by these 
pathogens can range from mild to very severe, with a mortality rate in some 
cases approaching 50%. 



A theory for the spread of t he Hantavirus was constructed a few years ago by 
two of the present authors ( Abramson and KenkreL 2002j ) and s hown to lead 



natur ally to spatiotemporal patterns such as the observed refugia (lYates et al 



2002 ) and the sporadic disapp earance and appearance of the epidemic ((Mills et al 
1999 : Parmenter et al. . 19991 ) . That theory will be referred to in the rest of 
the paper as AK. Seasonal, as well as extraordinary variations in demographic 
and environmental conditions are included in the AK model through spatio- 
temporal dependence of several parameter s, such as the car r ying c apacity. The 
emergence of traveling waves of infect ion ^Abramson et all M). the investi- 
gat ion of f luctuations (lAguirre et ah . 20021 ). extern al changes in environmen- 
tal ef f ects (iBal ard et alll2004l). and othe r features ( Kenkre , 20031 : Abramson , 
200l iKenkrel . l20oi " Bu ceta et all 120041 ) have also been studied. There are 
additional factors, not yet analyzed theoretically, which may be of importance 
in the dynamics of the virus. For example, human activity such as changes in 
agricultural practice may alter habitats and drive the rodent population into 
new habitats not previously occupied by them. While many theoretical issues 
regarding the AK development have been, and are being, explored quite in- 
tensely, the major problem of obtaining the values of the parameters inherent 
in the theory has remained neglected as a result of paucity of available data. 
Needless to say, the solution of this problem is crucial to the quantitative de- 
scription of the spread of the epidemic. The purpose of the present paper is 
such extraction of the essential parameters necessary in the description of the 
spread of the Hantavirus. 



The AK model (A bramson and KenkreL EooilAbramson et alll2003h is based 
on a fundamental set of biological features that characterize the transmission of 
Hantavirus among rodent populations, and involves five parameters: the birth 
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rate b of the rodents (mice), their death rate c, the environmental parameter 
K, the contagion rate a, and the diffusion coefficient D. In terms of these, the 
mice populations M s (susceptible) and Mj (infected) obey 



d -§ = HM, + Mi) - CM, - - aM,M, + DV>M„ (1) 

i)M, MAM S + Mi) 

All the parameters except X are considered to be independent of time t and 
space x. 

The measurement of the birth and death rates b, c presents no special chal- 
lenge. The environmental parameter K is sometimes measured from food and 
vegetation measurements, and also often obtained from aerial photographs of 
the landscape. Its variation in time and space can be well characterized al- 
though its absolute values are difficult to obtain. The encounter infection rate 
a is notoriously hard to measure from observations of individual mouse-mouse 

Table 1 

Summary of the mark-recapture data of Z. brevicauda. A total of 846 captures, 
corresponding to 411 different animals, were obtained. Of them, 188 were captured 
at least twice, and at most 10 times. The probability of recapture depends on sex 
and age group, as shown. It suggests an increase of the probability of recapture 
with age, independently supported by the probability of recapture as a function of 
weight (not shown). J: juveniles, SA: sub-adults, A: adults, F: females, M: males. 

J SA A Total 

Captured only once: 

F 7 13 98 118 



M 


8 


18 


79 


105 


Total 


15 


31 


177 


223 


Recaptured at least once: 










F 


2 


2 


71 


75 


M 





16 


97 


113 


Total 


2 


18 


168 


188 


Probability of recapture: 










F 


0.22 


0.13 


0.42 




M 


0.00 


0.47 


0.55 




Total 


0.13 


0.37 


0.49 
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interactions ( Botten et al. . 2002j ). At this stage of observational technique, we 



must assume that it is a floating parameter. 

The last of the five parameters, the diffusion constant D, is crucial to the AK 
description since the assumed mechanism for the spread of the epidemic is the 
diffusion (movement) of infected mice over the terrain followed by the trans- 
mission of infection to susceptible mice. In principle, it appears straightforward 
to measure D fr om records o f mice movement in a mark- recapture experiment, 



as was done by lOvaskainen (2004) in a recent study on the dispersal of but- 
terflies in a heterogeneous habitat. The investigation reported in the present 
paper began as an attempt to extract D directly in that manner. We will see 
that examination of the data has indeed allowed us to obtain values of D but 
also led us to a number of important conclusions about the transmission of 
Hantavirus among rodent populations, and has suggested substantial changes 
to be introduced into the theory of the spread of Hantavirus. 



2 The data set: Zygodontomys brevicauda in Panama 



At the beginning of 2000, human cases (more than 20 in the first cluster) 
of Hantavirus Pulmonary Syndrome (HPS) were recognized from the Azuero 
Peninsula, Panama. Following the outbreak, it was discovered that the pigmy 
rice rat Oligoryzomys fulvescens and the short-tailed cane mouse Zygodonto- 
mys brevicauda harbored two novel hantaviruses, Choclo virus (responsible 
for the HPS cases), and Calabazo virus (not known to cause human disease) 
respectively ( Vincent et al. . 2000f ). The data set we have selected in this study 



was obtained as the result of a mark-recapture observation performed in the 
Azuero Peninsula, in Tonosf (Los Santos), Panama, from June 27 to November 
20, 2003. The observation corresponds to the rainy season in the region, and 
following several years of relative draught. During these years, clinical cases of 
HPS were rare (after being more frequent in 1999-2001) and rodent densities 
were probably relatively low. 

Measurements were made on several species of rodents, of which we choose for 
the purposes of the present paper Z. brevicauda, host of Hantavirus Calabazo. 
This choice was made because Z. brevicauda was the most abundant species 
in the field study. A summary of the characteristics of the populations is 
presented in Table 1. The data relevant to the analysis of the movement consist 
of the position and time of capture of those mice that are captured at least 
twice, thus allowing for the calculation of their displacements from one location 
to another. The measurements were made with a square array of 7x 7 Sherman 
traps, separated 10 m each. Each measuring session lasted 3 days, with a time 
of recurrence to the same site of about one month. The number of trapping 
grids used in the study was 24, in 4 different sites. Each grid was set up across 
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Fig. 1. Displacements of the recaptured mice, projected along the x direction which 
is the direction into and out of the forest relative to the pastures in the landscape. 
Each point corresponds to a displacement of one of the 411 mice that were captured 
more than once. Some mice were captured more than twice, thus contributing with 
more than one point to this set. 



the edge between forest and pastures ([Suzan et all 120041 ) . For our analysis 
given below, all grids have been rotated such that the edge runs along the y 
direction. 



3 The movement of Z. brevicauda 



A few animals were captured a sufficient number of times (about 10), during 
a period of months, as to allow us to form a useful picture of the mouse walks. 
Unfortunately, there are only five of these mice. The limited number makes it 
impossible to carry out any statistical analysis of the properties of the walks. 

Despite this drawback, we were able to proceed with the analysis because 
the data set contains hundreds of recapture events, corresponding to different 
mice, each one providing us with a displacement at a certain time scale. These 
time scales are 1 and 2 days (if the recapture occurs during the same session, 
that lasts three days); about 1 month (if the recapture occurs at the next 
session); and about 2 months (at the second next session), 3 months, and 4 
months, respectively. At progressively longer time scales there are, certainly, 
less recapture events; but there are a sufficient number of them for sensible 
results to be obtained up to the 3-month scale. The data are shown in Fig. 1, 
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for the x-component of the 2-dimensional displacements of the mice. Each 
point in the graph represents a displacement Ax taking place during an in- 
terval At. Each displacement corresponds to the movement of a single mouse, 
but different displacements may or may not correspond to the same animal. 
We consider these displacements as our elementary events. Furthermore, we 
specifically assume them to be statistically independent. There are repeated 
events in the set, namely displacements with the same Ax and At, a fact not 
represented in Fig. 1. In the plot, it can be observed that, besides the 1- and 
2-day intervals, the data are scattered in "clouds" around 1, 2, 3 and 4 months. 
The reason for this is twofold. First, each session consists of 3 consecutive days 
of capture; therefore, the interval between two recaptures in different sessions 
is not a constant number of days. Second, variations due to the logistics of 
field work result in the time between sessions not being precisely 30 days. 

On each time scale, the set of available displacements is taken to represent 
a statistical ensemble, i.e., a population of ideal mice with certain statistical 
properties. Also, the displacements measured on each time scale correspond 
to a progressively coarser graining of the actual mice walks, containing an 
indeterminate (but presumably large) number of steps already on the 1-day 
scale. 

Using available data on each time scale, we construct mouse walks by ran- 
domly shuffling the displacements. These walks represent instances of possible 
walks, on that time scale, with exactly the same statistical properties as a 
hypothetical "representative mouse". We construct various such walks and 
perform ensemble averages to obtain the mean square displacement as a func- 
tion of time. For example, on the time scale of 1 day, 20 instances of the walk 
produce the result shown in Fig. 2. Diffusive behavior is inferred from the 
clearly linear rise of the mean square displacement. On longer time scales the 
same analysis produces smaller diffusion coefficients. However, application of 
the method is not reliable on these longer scales because it does not permit a 
correct evaluation of errors. 

The above considerations clarify the underlying diffusive nature of mice mo- 
tion. A direct analysis of the probability distribution function P(Ax) is pos- 
sible, also on each time scale. P(Ax) is a bell-shaped distribution, with a 
well-defined variance that can be used to characterize the evolution of the 
mean square displacement of the representative mouse as a function of time. 
The result of this analysis is in Fig. 3, where we show the mean square dis- 
placements in both the x and the y directions. Unlike in a simple diffusive 
process, the observed mean square displacement is observed to saturate to a 
value near 400 m 2 . The saturation value is different for (Ax 2 ) and for (Ay 2 ), 
indicating an anisotropy in the system. We have already mentioned above that 
each trapping grid contains an edge between forest and pastures, mainly along 
the y direction. This is surely responsible for the anisotropy. 
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Fig. 2. Square of the displacement as a function of time, on the time scale of 1 day, 
for 20 artificial "mouse walks" based on the recorded displacements. The heavy line 
is an ensemble average. 
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Fig. 3. Mean square displacement as a function of time, using the different time 
scales available from the data of Fig. 1. The two curves correspond to the two 
directions in space. 
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4 Movement in Confined Spaces: Rodent Home Ranges and Grid 
Sizes 



The saturation of the mean square displacement, clear from the observations 
discussed above, implies the existence of a spatial scale in the rodent system. 
One obvious possibility is that the saturation is a manifes tation of th e confined 
motion of the rodents, i.e., the existence of home ranges (Burt, 1943| ). In order 
to treat the problem quantitatively, we first present a simple calculation of the 
effect of confinement on the mean square displacement of a random walker. 



4-1 Effect of Confinement on (x 2 ) 



There are several ways one might model the effect of confinement on the mo- 
tion of rodents. One is to treat the motion as obeying, not a diffusion equation 



but a Fokker-Planck equation in an attractive potential ([Kuperman et al. 



2004). A simpler way, which we adopt here, is to consider the motion as oc- 
curring via simple diffusion but confined to a box whose size represents the 
home range. Restricting the analysis to 1-dimension for simplicity, we solve 
the diffusion equation for the probability per unit length P(x, t) of finding the 
mouse at position x and time t inside a bounded domain of length L. Taking 
into account the symmetry of the problem, the general solution can be writ- 
ten as P(x, t) = A + J] A\ cos [Xx] e~ x Dt . By imposing the condition that a 
mouse cannot escape from the home range (zero flux boundary condition), the 
allowed A coefficients can be calculated and the solution written as 



I 2 +°° 
P(x,t) = - + -J2A n cos 

L L n=l 



2mrx 



L2 



(3) 



wherein 



A, A 



,L/2 


'2mrx' 


I dx cos 


[ L \ 


l-L/2 





P{x,0). 



(4) 



With the initial condition that the mouse is within a length a centered around 
the origin of the box, specifically 



P(x,0) 



1/a \x\ < a/2 







|x| > a/2 



(5) 



the evolution at all times can be written as 



1 2 + °° sin (^) cos 



2n-KX \ „ „ 

— T — I (2n) /i iT /S Dt 

-e . 



an 



n=l 



(6) 
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The initial condition (5) reduces to a P(x, 0) = 8(x) in the limit a — > 0. The 
calculation of the mean square displacement (x 2 ) = f^^ 2 dxx 2 P(x,t) gives 



, , , +oc (-l) n+1 s infe) pjoW. 



07T "l n 



We see that (a; 2 ) ~ a 2 /12 + 2_Dt for t — > 0. We also see that (a; 2 ) saturates 
to L 2 /12 for i — > +oo. If P(rr, 0) = the mouse will move initially as if 

no home range existed. However, its diffusive motion will be limited in extent 
by the presence of the home range and (x 2 ) will eventually saturate. The 
1-dimensional calculation captures the essential features. The 2-dimensional 
extension, appropriate to mouse movement on the terrain, is straightforward 
to obtain because P(x,y,t) is given by the product of two functions of the 
form (6), one for the x-plane with D x and L x and the other for the y-plane with 
D y and L y . Allowing for the fact that differences in terrain could be reflected in 
the differences in D x and D y (given that the traps are laid out systematically 
relative to 'edges' in the landscape between forests and pastures), we give a 
usable expression for the average mean square displacement in 2-dimensions 
which can be used directly for the interpretation of the observations. 



<* 2 > + (y 2 ) 



L 2 X + L 2 L 1 +2?(-ir~ ! 
12 



l— 



X 



x { L. v sin I -j— J e Lx + /.,, >in y~J~ J ' 



Application of this analysis to the mouse data is straightforward. The short- 
time part of the mean square displacement gives the diffusion constant, av- 
eraged over the two directions, to be 200 ± 50 m 2 /d. The saturation value 
appears to imply that the home range L equals about 70 m. While these 
initial considerations suggest that the Panama data allow us to confirm the 
existence of rodent home ranges, as well as to measure their extent, careful 
observation introduces a note of caution: we notice that the derived value of 
the home range is of the order of the size of the measurement grid G = 60 m. 
Could sampling from a limited domain in space lead to saturation and mislead 
one into drawing incorrect conclusions about the home range? To answer this 
question we carry out the following analysis. 

4-2 Effect of Limited Spatial Observations on (x 2 ) 



Consider the motion of a random walker in unbounded space (no home ranges) 
but let the mean square displacement be calculated from observations in a 
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limited part of space of size G placed symmetrically around the origin for 
simplicity. We obviously have 



f' 2 dx x 2 P (x, t) 

( x ) = ^% — ( 9 ) 

Lg/2 dxP(x,t) 



where the probability P(x, t) is given by the propagator of the diffusion equa- 
tion in unbounded space: 



P( ^VS5f (10) 



Substitution gives 



(x 2 ) = 2Dt 



Ge isot 



V^erf(^) 



At short times (x 2 ) ~ 2Dt while at long times it saturates to the value G 2 /12. 
The 2-dimensional result is trivially obtained as a generalization of (11). 

Note that the behavior of a system without home ranges but with a finite 
window of observation is qualitatively different from that of a system with 
home ranges and an infinitely large window. The mean square displacement 
in the former (Eq. (11)) has quite a different time dependence from the latter 
(Eq. (7)), since the error function expressions differ considerably from expo- 
nentials. It is easy to see that in both cases the mean square displacement 
starts out at short times as 2Dt and saturates to a constant at long times. 
The saturation is to G 2 /12 in the first case and to L 2 /12 in the second. There 
is thus potential for confusion. One could mistakenly interpret what is actually 
the measurement grid size G to be the rodent home range L. 

Combination of the two elements discussed above, the home range effect and 
the grid size effect, is straightforward to analyze. Consider the situation in 
which a mouse is moving randomly inside a home range of width L but is 
observed only inside a region of width G, both concentric for simplicity. If the 
grid is larger than the home range (G > L), the evolution of the mean square 
displacement is exactly given by the previous result (7). In the other case, 
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Fig. 4. Mean square displacement in units of G 2 /12 as a function of time in units 
of G 2 /D for different values of the ratio G/L, i.e., the ratio of the measuring grid 
divided by the width of the home range. The dashed, solid, dash-dotted and dotted 
lines correspond, respectively, to G/L =1.43, 1, 0.7, 0. When G > L, (x 2 ) is seen 
to saturate to the lower value L 2 /12, whereas for all other (G < L) cases, it satu- 
rates to G 2 /12. The time dependence differs qualitatively for the various cases. The 
difference might be discernible — in principle — with more refined measurements. 

(G < L), the evolution is 



(x 2 ) 



91 

12 



■x 



+oo 



(2n)2-K z Dt 

e 



n=l 



oo 

. a r v Sin 

c + } 

= aw z L — / 



in(s22) gWriTrC) ( 2 ") 2 
^2 e £ 



n=l 



(12) 



wherein ( = G/L. Equation (12) reduces to Eq. (7) when ( — 1. It can 
be shown that (x 2 ) goes as a 2 /12 + 2Dt for t — > and saturates to L 2 /12 
for t — > +oo. In Fig. (4) we compare Eqs. (7), (11) and (12) for an initially 
localized condition P(x, 0) = 5(x). Despite the fact that all of the three curves 
(with G < L) show the same linear behavior at short times and reach the same 
saturation value, they differ considerably at intermediate times. A fourth curve 
(dashed, lowermost in the group) for which G > L is also shown in the same 
graph. 

An important step in the interpretation of saturation data is the verification 
of the sensitivity of the theoretical prediction to the initial position of the 
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Fig. 5. Mean square displacement in units of G 2 /6 as a function of the ratio of the 
home range size to the grid size, L/G. The curve is the result of a simulation of 
105 steps for each value of L/G. The arrow shows the average (in both directions 
x and y) mean square displacement observed in the measurement, from the data 
shown in Fig. 3, and the inferred value of the home range (L/G = 1.16). The grey 
lines show the error bounds. The dotted lines show the analytic result derived in 
the Appendix, valid for home ranges concentric with the grid. 

mice, as well as to the position of the home ranges with respect to the grid. 
In the Appendix we provide further details on the calculation of the mean 
square displacement, relevant to the dependence on the distribution of initial 
locations of the mice. To assess the combined effect, we have done a numerical 
experiment in which we measure simulated mice displacements. Each mouse is 
supposed to occupy uniformly its own home range, and all the home ranges are 
uniformly distributed in space. We limit the calculation of the mean square 
displacement to a window of size G, thus simulating the conditions of the 
measurement in the field. The result is shown in Fig. 5, where the normalized 
mean square displacement is plotted as a function of the normalized home 
range size L/G. The available data allow us to conclude that the home range 
of Z. brevicauda in Panama (averaged in the two directions) is 70 m, with 
asymmetric error bounds from 50 to 120 m. 



5 Discussion 

Extraction of quantitati ye information concerning the parameters in the A- 



bramson-Kenkre theory 
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was the initial task undertaken in the present investigation. This was prompted 
by the success of that theory in reproducing qualitatively observed features 
such as the sporadic disappearance of infection during time periods of low car- 
rying capacity, and the existence of refugia where infection persists and from 
where it propagates in the form of waves when conditions become favorable 
for such propagation. The key to the AK description of these processes was a 
transcritical bifurcation predicted by the theory and controlled by the carry- 
ing capacity. Crucial to quantitative application of the theory is the diffusion 
constant of the rodents, and the first aim of the present investigation was to 
extract this diffusion constant from movement observations collected for Zy- 
godontomys brevicauda in the peninsula de Azuero in Panama in a 5-month 
period in 2003. While l o gistic growth which is also part of the theory of ref 
(Ab ramson and Kenkrei . Effll is generally consi dered well established, both 
from field and laboratory studies ( Murray . 1993j ). equally strong justification 
for assuming diffusive transport for rodent movement is not available. The 
goal of extracting the diffusion constant has been met. The result particular 
to the species and location considered is D = 200 ± 50 m 2 per day. 



During the process of this extraction, we encountered a spatial scale in the 
rodent system which could be representative of the home range of the mice. 
Preliminary investigations showed that the observed spatial scale could be 
reflecting the size of the measuring grid instead of a characteristic of the ro- 
dent system. Since accuracy of interpretation will be considerably increased 
by having multiple values of the grid size, our analysis underscores the impor- 
tance for additional observations to be undertaken with grids of varying size. 
Motivated by this idea, we carried out further developments of the our theory 
which resulted in the simulation curve shown in Fig. 5. The mean value of the 
home range for the two directions can be read off from the plot to be L x =60 
m and L„=90 m. 



Because our prescription for extracting the diffusion constant relies on the 
initial time evolution of the mean square displacement rather than its satu- 
ration value, we have been successful in obtaining a usable estimate of the 
diffusion constant. However, even that evaluation suffers from the fact that 
reconstruction of mouse walks from observations is difficult because multiple 
measurements for individual mice are extremely rare. Therefore, we suggest 
strongly that systematic measurements of rodent movement be undertaken in 
other ways such as radio telemetry to confirm and sharpen the estimates of 
the diffusion constant we have obt ained in the present ana lysis. Elsewhere we 



will present our separate analysis (lAbramson et all 120041) o f mice movement 



data taken from web measurements ( Parmenter et all l2003[ ). 



The present analysis has pointed to the importance of the concept of home 
ranges, which have been known and discussed in many places in the literature 
earlier (Burt, 19431 : Anderson . 19821 : Fordl 1979f ). They would generally be re- 
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fleeted in animal movement measurements in the way we have detailed in our 
calculations in the present paper and could additionally be of crucial impor- 
tance in the theory of the spread of epidemics. While they do not appear in the 
AK description (the tacit assumption in that formalism being that they are 



large r than other lengths of interest), we have now developed ((Kenkre et al 



2004) a series of models to treat them explicitly. The basic idea in our new 



theory is to consider the dynamics of two types of mice, stationary and itin- 
erant (and susceptible and infected in each category). The stationary mice 
are the adults that move within their home ranges and do not stray far from 
the burrow. The itinerant mice are the subadults that must leave to find their 
own home ranges. Our studies employ a combination of nonlinear analysis and 
simulations. The results will be presented elsewhere. 
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A Dependence of (x 2 ) on distribution of initial location of rodents 
within the grid 



The mean square displacement obviously depends on the initial location of 
the mice, x . Let us consider the case of infinite home range, but finite G, and 
determine the effect on (11) of the distribution V(x ) of initial locations of the 
mice. We thus evaluate ((x — xq) 2 ) and calculate its average, denoted by an 
additional (•••), with respect to the distribution V(xq) of the initial positions 
through (((x — xq) 2 )) = dx V(xo)((x — x ) 2 ). Proceeding as in Eq. (11), 

we have 
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This reduces to Eq. (11) when V(x ) = 5(x ). The evolution at short times is 
found to be independent of V(x ): 



hm(((x-x ) 2 ))~2Dt, 



(A.2) 



while the saturation turns out to depend on V(x ) and is given by 
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Notice that for a uniform initial distribution V(xq) = 1/G, the integral equals 
G 2 /12, and the saturation value is G 2 /6. This value is twice that obtained for 
the case of initial placement of the mice at the center of the grid. We see from 
(A.3) that its right hand side lies always between G 2 /12 and G 2 /3 depending 
on the initial distribution. 



The time-dependent evolution of (((x — x ) 2 )) can also be determined when 
the home range is not infinite. If Xq is not the center of the home range, P(x, t) 
is no longer symmetric with respect to the origin, and the series solution for 
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P(x,t) now contains sine (in addition to cosine) functions: 
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Here the initial probability is 

Px o (x,0) = < 



1/ct \x — xq\ < a/2 
\x — x \ > a/2 



(A.5) 



We continue to consider the home range and grid window to be concentric for 
simplicity. The calculation follows the steps shown in earlier cases and yields, 
for ({(x — x ) 2 )), the expression: 
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The time dependence of Eq. (A.6) at short and long times can be shown to 
tend to that given by Eq. (A. 2) and (A. 3), respectively, for the case of free 
diffusion. In the case L < G, the result is simply Eq. (A.6) with the parameter 
C set equal to 1 and G set equal to L. In such a case the saturation value is 
given by an expression similar to Eq. (A. 3) but with L replacing G everywhere. 



18 



